Sterile neutrino signals from supernovae 
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We investigate the effects of a mixing of active and sterile neutrinos on the ratios of supernova 
electron neutrino flux [Fe) and antineutrino flux [Fe) to the total flux of the other neutrino and 
antineutrino flavours {Fa). We assume that the heaviest (in the normal hierarchy) Standard Model 
neutrino mixes with a sterile neutrino resulting in a pair of mass eigenstates with a small mass 
gap. Using the density matrix formalism we solve numerically the the evolution of neutrino states 
in the envelope of a supernova and determine the flux ratios F^/Fa and F^/Fa as a function of the 
active-sterile mixing angle and for the experimentally allowed range of the standard active-active 
mixing angle 613. 

PACS numbers: 14.60 Pq, 13.15 +g, 95.85 Ry 
Keywords: sterile neutrino, neutrino mixing, supernova 



Oh' 
I 

Oh' 
CU 



(N 
> 

cn 
m 
rn 

00 

o 
o 



X 



Introduction. The neutrino oscillation phenomena ob- 
served so far can be explained in terms of three active 
neutrinos V(.,v^ and (for reviews, see [T*], [5]). The 
only exception to this has been, until recently, the muon 
and electron neutrino and antineutrino oscillation claim 
from LSND experiment [sl, whose consistency with the 
other data would have required the existence of a fourth 
neutrino type, a sterile neutrino. The recent results from 
MiniBooNE experiment [1] seem to refute this claim, 
however. Hence, as far as the interpretation of the ex- 
isting experimental data is concerned there seems to be 
neither need nor too much space for sterile neutrinos Q . 
(Notice, however, Ref. (H). 

Nevertheless, it is quite possible that sterile neutri- 
nos do exist anyway. Many electroweak theories beyond 
the Standard Model predict the existence of right-handed 
neutrinos, putting neutrinos in this respect on the same 
footing as the other basic particles, charged leptons and 
quarks. If the right-handed neutrinos exist they would 
be sterile under the SU(2)xU(l) interactions of the Stan- 
dard Model (SM). Many Grand Unified Theories, e.g. the 
SO(IO) model, include sterile neutrinos, too. 

Obviously, sterile neutrinos would have escaped de- 
tection if their mixing with the active neutrinos is very 
small or iron-existing. Another reason for their elusive- 
ness could be that they are closely degenerate in mass 
with the active neutrinos and therefore kinematically in- 
distinguishable from them. The two close-mass eigen- 
states resulting from the mixing of an active neutrino 
and a sterile neutrino could have appeared in experiments 
performed so far as a single neutrino with the normal SM 
properties, not as two separate states with interactions 
deviating from the SM. 

We have shown earlier that degenerate pairs of neutri- 
nos resulting in mixing of an active and a sterile neutrino, 
while not separable in laboratory experiments or in the 
oscillation phenomena studied so far, might reveal them- 
selves by affecting the flux rations of the active neutri- 
nos Ve, I'll and Vt measured at the Earth and originating 



in distant astrophysical objects like active galactic nuclei 
(AGN) and supernovae [8[ . Neutrino fluxes from these 
sources are sensitive to smaller squared mass differences 
Am^ than those detected previously. 

In this paper we will return to the subject and study in 
more detail the effects of sterile neutrinos on the flux of 
electroir neutrinos {F^) and electron antineutrinos (i^g) 
from supernovae. We have developed a numerical code 
that computes the evolution of neutrino energy states in a 
medium of varying density and determines the flux ratios 
Fe/Fa and F^/Fa, where Fa is the total flux of muon and 
tau neutrinos and antineutrinos. Fa = Ffj_ + F11 + Fr + Ff. 
We find out that even quite a small mixing of a sterile and 
an active neutrino might make the ratios F^/Fa and in 
particular Fg /Fa to deviate substantially from the values 
they obtain in the case of three active neutrinos. 

The scheme. We will work in the following active-sterile 
mixing scheme Q- We define 



(* = 1,2,3), 



(1) 



where U is the standard unitary 3x3 matrix parame- 
terized by three rotation angles ^12,^23 and ^13 and 
{£ = e,fi,T) are the ordinary active neutrinos. We then 
assume that each state is accompanied by a sterile 
state I'si, and that these two states mix (to be called as 
an zi'-mixing) according to 



Ui = cos (f)i£'i - sin (l)th'si, 
v[ = sin (jj^iyi + cos 4)^1' si, 



(2) 



and form two mass eigenstates with masses and m^, 
respectively. In the limit where the mixing angles 4'i or 
the squared mass differences Am|j, = mf — m/^ , or both, 
are very small, the rotation angles 6ij approach the ro- 
tation angles 6*12, 6*23 and 6*13 of the standard three ac- 
tive neutrino case and I'i's the three ordinary mass eigen- 
states. 

In this work we will use the following experimental best 
fit values for the Standard Model mixing angles 9i2 and 
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6*23 and the squared mass differences [2|: 

9x2 - 33.7°, 

02i - 43.3°, 

Ato?2 = 7.9•fO~^eV^ 

Ato^3 = 2.6-fO-^eV^ (3) 

For the mixing angle ^13 we use the value ^13 = 5.2°, 
which is the experimental upper limit at la level ( at 3(7 
level ^13 < ff.5° ) 0. 

Solar constraints. Before moving to supernova physics 
we have to study the restrictions the solar neutrino data 
sets on the active-sterile neutrino mixing we are consid- 
ering. The SNO solar neutrino experiment [3] has es- 
tablished boundaries for the relative neutrino fluxes of 
the active neutrinos arriving from the Sun Sterile 
neutrino mixing would affect these ratios, and by re- 
quiring that the effects are within the data uncertain- 
ties one can constrain the values of active-sterile mixing 
parameters 4'i and Am^j/ . Our numerical study shows 
that solar neutrino data restrict the 11'- and 22'-mixings 
strongly leaving room only for mixings with very small 
Am^ or very small mixing angle, which is in good ac- 
cordance with earlier work done by other authors (e.g. 
[lOj). The restrictions arise mainly from the low energy 
part of the flux, where the vacuum oscillation dominates 
strongly over the matter effects. It turns out that given 
these tight constraints, the 11'- and 22'-mixings would 
not have noticeable effect on supernova neutrino fluxes. 
For the 33'-mixing, in contrast, there is practically no re- 
strictions from the solar neutrino data. This is due to the 
almost negligible Vg, proportion in the and ^ states. 

This leaves a possibility for large Amgg/ and ^3 , which 
allows the active-sterile mixing to have a detectable ef- 
fect on supernova neutrino fluxes as we will show in the 
following. From now on we will assume that sterile states 
Vsi and effectively decouple and that only the sterile 
state affects the behaviour of active neutrinos. 

Numerical study. We determine the evolution of the 
states of neutrinos traversing the supernova envelope by 
solving numerically the density matrix equation 

^=*[p,i?], (4) 
where the Hamiltonian H is given by 

H = ^Udiagimlml ml m'i)U^ + diag(y, 0, 0, V/2), 

(5) 

where V = ^/2GfNi,, and is the number density of 
electrons {llj , assuming that the neutron density is equal 
to the electron density. The evolution of antineutrino 
states is obtained similarly, except that the matter term 
V in the Hamiltonian has the opposite sign. 

The flavor decomposition of the neutrino flux as a func- 
tion of time is given by the diagonal elements pu of the 
density matrix once the initial fluxes are given. 



We will use the typical density profile of a progenitor 
star as given in 12]. This yields the potential (r ~ t) 

F(r) = 1.5- 10^2 eV(m/r)^ (6) 

Following [1^ we assume that the progenitor star extends 
to i? = 2.8 • 10^° m. That is, we take r to vary from the 
neutrino sphere at about 10^ m to R. 

In Fig. 1 the energy states Vi {i = 1,2,3) and z/3 are 
plotted as a function of the electron density inside the 
star, starting from the surface, where the density of mat- 
ter and subsequently the potential V become zero. At 
high density the potential V dominates the Hamiltonian 
([5]) and the energy eigenstates are close to the flavour 
states vi and ly^ {£ = e, /i, t, s). At zero density the states 
equal to the mass eigenstates 



H 




FIG. 1: Neutrino eigenstates in matter in the presence of 
33'-mixing. The L-resonance is always adiabatic, whereas be- 
havior in the H- and 33 '-resonances depend on the value of 
the corresponding mixing angles S13 and Lp. 

We assume the initial neutrino and antineutrino fluxes 
to be as given in [Tij. The authors of that reference 
have developed a model for supernova explosion and 
based on this model they give the fluxes of I'e, i^e and 
j/^ ~ i/^ ~ j/^ ~ = as a function of neutrino en- 
ergy. More specifically, we will use in our analysis the 
ratios : Fg : = A : 2.3 : 1.4, which correspond to 
the fluxes at the energy 10 MeV where the intensity of 
the neutrino radiation is largest. Of course, the fluxes, 
-F!j? in particular, will be hard to measure within a sin- 
gle energy bin [16]. Nevertheless, the flux ratios at the 
peak energy 10 MeV turn out to be very close to the in- 
tegrated fluxes, with just slightly overweighting the frac- 
tion of electron neutrinos, so that they can be reliably 
used as input fluxes. Sterile neutrinos are of course not 
produced in supernovae but they are created, as a result 
of their assumed mixing with the standard model neutri- 
nos, during the traverse of neutrinos the medium of the 
star. 

Results. Let us first assume that sterile neutrinos are 
absent and study the system with the three active neu- 
trinos and their antineutrinos only. By using these as 
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initial conditions we solved numerically the density ma- 
trix equation and determined the the flux ratios F^/Fa 
and Fg/ Fa- Allowing the mixing angle to vary in the 
range < 6*13 < 5, 2° and fixing the other mixing pa- 
rameters to their best-fit values ^ the flux ratios in the 
case of three active neutrinos are found to vary within 
the ranges 



0.17 < FjFa <0.29, 
0.24 < F^/Fa <0.26. 



(7) 



The sensitivity of these ratios on the mixing angle ^13 
is due to the fact that the adiabaticity of the so called 
i?-resonance at the crossing of the energy states V2 and 
1^3 (see Fig. [T]) depends on the value of this angle. When 
^13 = 0, the behavior of the system is fully non-adiabatic, 
meaning that the flux of electron neutrinos, which is the 
largest one initially, ends to the energy state 1^2, which 



contains less of Va (= v^j_ 



Vr) than the energy 



state V'i. This corresponds to the largest values of the 
ratios F^/Fa and Fg/Fa. The smallest values of these 
ratios are obtained when ^13 attains its largest exper- 
imentally allowed value. In this case the if-resonance 
is adiabatic and the initial electron neutrino flux ends 
into the state z/3, which has a very small i>e component, 
|(i/e|i^3)|' = sin2 0i3 = 0.008. 

We now add into the scheme a sterile neutrino. We 
assume that the sterile neutrino mixes with the active 
neutrinos in such a way that the states vi and z/2 remain 
intact, that is, their flavor decomposition is the same as 
in the Standard Model and consist only of active flavors, 
while the third state 1^3 becomes a superposition of the 
active neutrinos and the sterile neutrino according to 
The state 1/3 is assumed to be close in mass to the new 
heavier state v'^ orthogonal to it. In our analysis we will 



assume m 



'2 



= 10^^ eV^. Increasing ATO33, from 
this value would not change our results. For lower val- 
ues of Am33, , the effect of active-sterile mixing decreases 



and will be washed out at Ar 



10-"eV% where our 



results would coincide with the SM result. As 1^3 is heav- 
ier than 1/3 , the energy levels corresponding to these two 
states cross in the antineutrino sector, not in the neutrino 
sector (see Fig. [1]). 

In Fig. [2] we present the flux ratios F^/Fa and Fg/Fa 
as a function of the active-sterile mixing angle (f> for the 
experimentally allowed range of the angle 9i3, < ^13 < 
5,2°. In what follows we will explain the behavior of the 
ratios qualitatively. 

The both ratios are larger for ^13 — 0° than for 613 ~ 
5.2°. This feature is easy to understand in the following 
way. When ^13 = 0, the i/-resonance in the neutrino 
sector is non-adiabatic and none of the electron neutrinos 
created in the supernova core ends up to the state V3 and 
the state ^3 will be populated solely by the less intense 
i^r component of the initial flux. The original electron 
neutrinos end up to the state 1^2 whose component is 
much larger than that of v^, where the original electron 
neutrinos will end to in the case of ^13 = 5.2. As to the 
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FIG. 2: Ratio of (a) (b) i/e and nonelectronic active fluxes 



as a function of sterile mixing angle (p. Results correspond- 
ing to the minimum and maximum values of 613 allowed by 
experiment are shown. The region denoted by SM is the al- 
lowed range in the SM case of three active neutrinos. Here 
the neutrino energy is taken as i? = 10 MeV. 



ratio Fe/Fa, the flux Fg does not depend on the value of 
6*13, but Fa is larger in the case of ^13 = 5.2° than in the 
case of 6*13 = 0° as then the original electron neutrino 
flux ends to the state , where there are more muon and 
tau neutrinos than in the state V2. 

The 33 '-mixing affects the ratios F^/Fa and Fg/Fa in 
two ways. Firstly, the mixing makes the flavor decompo- 
sition of the energy eigenstate 1^3 to change from its SM 
form as it replaces the active flavors partly by the sterile 
flavor. As there is only a tiny electron neutrino (antineu- 
trino) component in 1/3 (P3), this would affect mainly 
the flux Fa and makes the value of the ratios F^/Fa and 
Fe/Fa larger than in the Standard Model case. The effect 
of course increases with increasing active-sterile mixing 
angle (p. This explains the growth of the ratios when (p 
shifts from zero to about 12° . 

The second effect of the active-sterile mixings is related 
to the possible level crossing in the resonance transition 
region of the D3 and v'^ antineutrino states. For a small 
mixing angle cj), the 33 '-transition is non-adiabatic and 
the level crossing takes place. That is, all Vt's created in 
the supernova will end up to the state V3 , which has only 
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a tiny sterile content, proportional to sin^ (f>. 

When (f) > 12° the 33'-resonance is fully adiabatic. In 
this case the original Vt flux will go into the state v'^, 
which has a sterile fraction proportional to cos^ (f>. When 
(j) increases, the and z/^ components of v'^ will become 
larger, which will make Fa to increase. What will hap- 
pen to the ratios depends, however, on the adiabaticity 
of the H-resonance. When 613 = the H-resonance is 
non-adiabatic, the occupation of the state 1/3 originates 
from the original Vr flux. When the angle (p increases, 
the and components of will become smaller and 
the sterile part will become larger. This will make Fa to 
decrease. As the original flux of tau neutrinos and an- 
tineutrinos are assumed to be equal, the net effect will 
be that Fa remains practically constant for (j> > 12°. As 
the electron neutrino and antineutrino fluxes are quite 
insensitive to the active-sterile mixing, this implies that 
the values of the ratios F^/Fa and F^/Fa saturate to the 
constant values for (f) > 12° and 6113 — 0. According 
to Fig. [21 the saturation values are Fg/Fa — 0.35 and 
Fg/Fa — 0.32, to be compared with the Standard Model 
values 0.29 and 0.26, respectively. 

When 9i3 = 5.2°, the H-resonance is adiabatic. Neu- 
trinos created originally as Ve will in this case end up 
to the state 1^3. In contrast with the previous case, now 
the ratios F^/Fa and Fg/Fa will not saturate after the 
active-sterile mixing has achieved the value (p = 12° but 
will keep increasing with (j). The difference to the pre- 
vious case is that now the state 1^3 will be more occu- 
pied than the state 1/3 as the original electron neutrino 
flux is assumed to be larger than the tau antineutrino 
flux. For the maximal 33'-mixing, i.e. for tp = 45°, one 
has Fe/Fa = 0.24 and Fg/Fa = 0.35, to be compared 
with the SM values 0.17 and 0.24, respectively. Note 
that Fg/Fa becomes insensitive to the angle 6113 when 
the active-sterile mixing angle has the value about 35°. 

Discussion and summary. As one can infer from the Fig. 
[21 the ratio Fg/Fa would give more information on the 
active-sterile mixing than F^/Fa taking the uncertainty 
of the value of the mixing angle ^13 into account. Given 
the value of Fg/Fa from measurements, the angle 4> will 
be limited in general into quite a narrow range of values. 
Also, for 0^5° the Fg/Fa would deviate from the SM 



prediction for all allowed values of ^13, in contrast with 

Fe/Fa. 

As for the uncertainties in the initial neutrino fluxes, 
the differences between the SM case and sterile case arise 
only due to the 33'-mixing in the antineutrino part. So 
only the initial flux of Dr will affect the ratios by possibly 
changing the value of Fa. 

It is worthwhile to make a remark concerning the effect 
of the Earth on the fluxes. As there is only a tiny elec- 
tron neutrino component in the state 1^3, the effects of 
neutrino-matter interactions in the Earth will not differ 
much from those in the SM case. The active-sterile mix- 
ing will only slightly change the fluxes of non-electron 
flavour neutrinos measured at the opposite side of the 
Earth. The most visible effect would be a distortion of 
the vi, {I e,e) energy spectra, which could possibl y b e 
measured in future supernova neutrino experiments [l5| . 

In summary, we have considered the situation where 
neutrinos appear as three pairs of states, where each 
pair orginates in a mixing of an active neutrino and a 
sterile neutrino and its components are so close in mass 
that they are not distinguished as different particles in 
the present experiments. We have shown that this kind 
of structure might manifest itself through the fractions 
of electron neutrinos and antineutrinos in the supernova 
neutrino flux, i.e. through the ratios F^/Fa and Fg/Fa. 
Apart from the active-sterile mixing angle 0, the devia- 
tions of these ratios from their Standard Model values are 
sensitive also to the mixing angle ^13 for which an upper 
limit only is known at the moment. Only an active-sterile 
mixing associated with the heaviest mass state 1^3 can 
have measurable effects, since the solar neutrino data re- 
strict the mixing parameters in the case of the two other 
states too small to lead any observable changes in the 
supernovae neutrino fluxes. 
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